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Summary 

The studies oarried out during the period from 15 August 1984 through 
30 November 1984 are dealt with in this Quarterly Report. Such studies are 
related to the following areas of investigation: 

a. The dynamic response of tho TSS during the entire station-koeping 
phase for the first electrodynamic mission. From the simulations that we 
ran useful information such as the out-of-plane swing amplitude and the 
tether's bowing have been derived. 

b. The dynamics of the 6lack tether. Our in-house high-resolution 
computer code, SLACK2, has been improved (additional work will be carried 
out in the future) both in capabilities and computational speed. A 
convincing test case has been run together with some fairly long 
simulations of a severed tether with variable longitudinal oscillation 
damping. 


c. Safety hazard studies related to tether breakage or plasms 
contactor failure. Preliminary values of tho potential difference after 
the failure and of the drop of the electric field along the tether axis 
have beeu computed. 

d. The up-date of the satellite rotational dynamics model. Such 
up-date has been initiated but results are not yet available (they will be 
shown in the next Quarterly Report). 


Figure Captions 


Figure A, B, C 


Figure D, E a F 


Figure G 
Figure H, I, J 


Figure K-l 


Figure K-2 


Figure K-3 


Figure K-4 


In-Plane vs. Time (MKS Units) 

Out-of-Plane vs. Time (MKS Units) 

Current vs. Timo (HKS Units) 

(The current must be doubled to bo 
consistent with the dynamic response.) 

1st Case of TSS First Eleotrodynamio Mission. 
Station-Keeping Phase. No Damping. 

In-Plane vs. Time (HKS Units) 

Out-of-Plane vs. Timo (MKS Units) 

Current vs. Time (MKS Units) 

(The current must be doubled to be 
consistent with the dynamic response.) 

2nd Case of TSS First Electrodynamic Mission. 
Station-Keeping Phase. With Damping. 

Geometry for the Tether Bowing Analysis. 

In-Plane Tether Bowing vs. Time (MKS Units) 
Out-of-Plane Tether Bowing vs. Time (HKS Units) 

Current vs. Time (MKS Units) 

First Electrodynamic Mission. Average Current in 
the Tether- 0.85 Amp. The Pendular Libration has 
been removed from the plots to better show the tether 
deflection wrt. the Orbiter/Satellite line of sight. 

The radial components vs. time for (a) a SKYHOOK rnn 
and (b) an equivalent SLACK2 run. See discussion in 
text. 

Successive configuration plots for the case considered 
in text (the SKYHOOK results are shown). 

The configuration plots for SKYHOOK and SLACK2 
simulation of same case (as discussed in text) are 
shown sido by side for several values of T. SKYHOOK 
results are on the left, inverted to make the in-plane 
components agree. 

Results from a SLACK2 run with a long tether remnant 
after break (20 kilometers remaining from 100 km 
tether; 30 segments iu simulation). This demonstrates 
that a small randomization of the initial conditions 
allows successful simulation in long remnant cases. 
Configurations are shown at 50 second intervals for 
2500 seconds. 


Figure K-S 


Figure K-6 


Figure L 
Figure H 


Figure N 


SLACK2 results for a sequence of runs in which only the 
tether damping varies. Otherwise, the cases arc all 
the same: 0.2 Jem remaining from a 20 Jem upward 

deployed tether; 35 segment discretization; 600 
seoond simulation with output at 10 second intervals. 
Damping factors of 0%, 0.5%, 1% and 2% are used, 
ac discussed in the text. 

The final case of Figure K-5 (2% damping) is followed 
for 3000 seoonds, showing configurations at 20 second 
intervals. 

Simplified equivalent circuit of electrodynamic tethor. 

In-vacuo, elongated prolate spheroid model of long 
orbiting tether. 

Geometry for electric field computation. 
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1.0 Introduction 


This is the first quarterly report submitted by SAO under contract 
NASB-36160, "Tho Investigation of Tethered Satellite System Dynamics," Dr. 
Enrico Lorenzini, PI, and covers the period from 15 August 1984 through 30 
November 1984. 


2.0 Technical Activity During Reporting Period and Program Status 

2.1 TSS Tether Dynamics and Out-of-Plane Swing Amplitude During 

the Station-Keeping Phase of the First Electrodynamic Mission 

2.1.1 Satellite Libration Dynamics 

In order to determine the out-of-plane swing angle that can be 
expected in an electrodynamic mission, an integration lasting 80,000 
seoonds has bean done for a case similar to what can be expected during the 
first flight. In the simulation a 100 metric ton Shuttle is at an altitude 

of 295 km with a 550 kg subsatellite deployed upward on a 20 km tether. 

Both the Shuttle and the subsatellite have plasma contactors and the 
resistance of the wire is 4000 ohms {value provided by MMA) . Since the 
integration is a long one only two mass points have been used in the 
simulation. In this model the eleotrodynamic force is applied at the end 
of the wire instead of being distributed along the wire so that the torque 
causing librations of the tether is too large by a factor of two. The 

results can be scaled by multiplying the current by a factor of two to give 

physically meaningful results (further details on this point will be given 
later in this report). 


The only properties of the wire modelled in the simulation are its 
stiffness and damping. The parameters for a 2 mm kevlur wire (the axial 
stiffness is supposed to be provided by the kevlar only) have been used 
which gives a stiffness of 10998 dyncs/om. A oriticol damping coefficient 
of 155122 dynes per cm/sco is used to suppress longitudinal vibrations that 
would cause slow numerical integration. Since the wire is stiff and no 
reel control algorithm is usod. thore is no damping of pendular 
oscillations. In the first run, the plasma contactors were modelled by 
using large balloon radii for the Shuttle and tho subsatellite. This has 
tho effect of grounding both ends of the wire to the plasma. Tho 
integration failed after about 15,000 seconds and the integration started 
taking very smnll integration steps. Electrodynamio integrations are 
difficult because of the very fast time constants of the electrodynamio 
variables. The model for charge collection by the balloon is non - linear 
and has discontinuities. A large balloon radius also makes the integration 
more critical. In order to make the integration easier, the balloon model 
was replaced by a one ohm resistor between the electrodes and the plasmas. 
This eliminates the non-linearity and discontinuities and accomplishes the 
objective of grounding the electrodes to the plasma. 

The integration was started with the system on the x-axis at the 
ascending node headed Easterly (toward the +y axis) in a 28° inclination 
orbit. The epoch of the orbit was chosen so that the north magnetic pole 

is inclined toward the -y axis. In this orientation the orbit has an 

o 

inclination of about 16.5 with respect to the magnetic equator. During the 
integration the poles rotates Easterly through almost one complete 
revolution. This rotation of the pole is the major factor affecting the 
electrodynamic force on the wire. 
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During the first half orhit (see Figure A, D and C), the in-plane 
libratiou angle reaches .6 Jan (1.7 deg) and the out-of-plane amplitude 
reoohe6 .15 bn (.4 deg). During the first 25.000 seoonds the out-of-plano 
amplitude builds up to an amplitude of about .9 bn (2.6 degrees) while the 
in-plane amplitude stays about tho same. After 25.000 sec there is no 
significant change in the ont-of-plane amplitude. Tho wire current 
averages about .85 amps over the rnn. This is equivalent to the torque of 
a current of 2 z .85 » 1.70 amps when the electrodynamic force is 
distributed along the wire instead of being applied to the end mass. By 
dividing the libration amplitude by a factor of two to suppress the 
overestimated torque, and by comparing it to the atmospheric mission case, 
it turns out that the out-of-plane angular momentum at tho beginning of 
retrieval, for the first electrodynamic mission, is comparable to the one 
evaluated for the atmospheric mission (100 km tether and 0.2 dog 
out-of-plane angle). The build-up in the out-of-plane oscillation 
amplitude is tho result of resonance betwoen the electrodynamic driving 
force which has a resonant (small) component and the natural period of the 
out-of-plane oscillation which is ot half the orbital period. The system 
can go off resonance because of a change in the period of either the 
driving force or the natural out-of-plane oscillation. The natural period 
lengthens with increasing amplitude, being exactly half the orbital period 
only for small amplitudes. The amplitude is still fairly small in this 
run, but the dependence of the period on amplitude has not been studied. 
Since the magnetic field is rotating with tho earth, the ohange in the pole 
position may be a more likely cause of the system going off resonanoe. 
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Table 1 below lists the moan value and the r,m. s. deviation of various 
quontlties for the run, along with the rescaled value of the current that 
actually produces this dynamio response. 


Table 1 



Average 

R.M.S. 

Minimum 

Maximum 


Value 

Value 

Value 

Value 

In- Plane 

Di spl acement 

186 m 

242 m 

-315 » 

+718 m 

In- Plane Angle 

0.53 dog 

0.69 deg 

-0.9 deg 

2.06 deg 

Out-of-Plane 

Displacement 

2.7 m 

495 m 

-897 m 

+930 m 

Out-of-Plano 

Angle 

0.01 deg 

1.42 dog 

-2.57 deg 

2,67 deg 

Rescaled 

Current 

+ 

-1.70 A 

.34 A 

-.76 A 

-2.32 A 

VxD. Z 

-3403 v 

696 V 

-1532 V 

-4640 V 


Table 1. Summary of the dynamio response during the 

entire station-keeping phase (cleotrodynamlo 
mission) of a TSS system without damping, 

Tho out-of-plane oscillation could be theoretically reduced in a 
couple of different ways. Probably the most effioient would be the use of 
electrodynamic damping since it is a conducting tether. This technique is 
discussed in detail in the report "Investigation of Electrodynamic 
Stabilization and Control of Long Orbiting Tethers," G. Colombo, Interim 
Report, NASA Contract NAS8-33691, March 1981. The electrodynamic damping 
can be toned to the out-of-plane oscillation while ignoring the effect on 
the in-plane oscillation which is easily handled by reel control techniques 
(as later shown). The eleotrodynamio damping can be applied during the 
retrieval along with a reel control technique. 


V 
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To Investigate the effect t of damping by reel oontrol, a similar, long 
duration simulation was run. Tho tether tension control law makes the 
tether respond an a highly visco-elastic medium. For this reason tho 
control law was simulated by assuming an artificial tether stiffness and 
damping coefficient that give the samo longitudinal dynamic response 
provided by the reeling mechanism. 

In order to maximize the ont-of-plane damping effectiveness the 
longitudinal oscillation frequency was fixed at twice the out-of-plane 
swing frequency Zonr times the orbital frequency) and the damping 
coefficient was chosen as low as £« 0.2. This low value of tho damping 
coefficient is derived from an in-house preliminary study oarried out at 
SAO. Tho lower the damping coefficient the larger the energy transfer 
between the pendular libration and the longitudinal oscillation. However a 
larger tether length variation corresponds to a smaller damping coefficient 
so that, the attainable values are lowerly limited, K “0.2 appears to be a 
good compromise in order to avoid excessive tether lengthening d<'e to the 
actual libration amplitude during station-keeping. In this simulation the 
overestimated effect (by a faotor of two) of the electrodynamic torque was 
compensated by halving the current in the tether. The results obtained 
represent therefore the worst case dynamic response of the system with the 
tether resistance presently planned by MMA (tho actually achievable ourrent 
in twice as ouch the one shown in Figure F) . In-plane response, 
out-of-plane responso and current are plotted in Figure D, E and F 
respectively. As expected the values of the oscillation amplitudes are 
very close to one half those obtained in the other simulation (except the 
in-plane that is affectod by the damping). This moans that the dynamic 
response (at least for small osoillation) is linearly dependent on the 
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ourrent; the model adopted for tho clootrodynamio torque is therefore 
fully acceptable. Results from the simulation show that the out-of-plano 
oscillation is almost unappreoiably affected by the damping (as suspected) . 
The out-of-plane damping depends, by the way, on the square of the 
out-of-plane oscillation amplitude and it is therefore very poor for 
amplitude around one degree. The in-plane, on the contrary, bonefits from 
the damping (even if the tether control is tuned to the out-of-plane) so 
that tho in-plane transient oscillation caused by the current activation 
(step funotion) is reduced, after 8 hrs., to its steady state value. A 
tether control tuned to the in-plane osoillation (as presently done by MHA) 
is oven more effective to damping out the in-plane oscillation. It is 
therefore more convenient to handle tho in-plane libration by tetbor 
oontrol and to damp (if deemed neoessnry) the out-of-plane libration by 
current control. In Table 2 the most important results of this simulation 
are summarized. In this table the current is rescaled (multiplied by a 
factor of 2) to the value that ootually produces the dynamic response shown 
in tho upper linos. 



Average 

Value 

Table 2 

R.H.S. 

Value 

Minimum 

Value 

Maximum 

Value 

In- PI ane 

Displacement 

88 m 

110 in 

-280 m 

477 m 

In- PI ane Angle 

0.25 deg 

0.32 deg 

-0.8 deg 

1.37 deg 

Out-of-Pl ane 

Displacement 

0.9 m 

239 m 

-433 m 

451 a 

Out-of-Plane 

Angle 

2.6°xl0" 3 

0.68° 

-1.24° 

1.29° 

Rescaled Current 

-0.84 A 

0.16 A 

-0.38 A 

-1.14 A 


Table 2. Summary of the dynamic response during the 

entire station-keeping phase (electrodynamic 
mission) of a TSS system with positive damping. 
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Note that the current adopted is the maximum attainable with the present 
tether reslatanoe and plasma contactors at both tether ends (satellite and 
Orbiter). The current in the tethor is resittively limited and therefore 
depends on tho interaction between the tether and the magnetic field only 
(no activo control). The current is switched on abrnptly (step function) 
and the dynamic response reprosonts, therefore, a worst case condition. 
Since tho out-of-plane dynamics is actually undamped the overshoot of tho 
system after the step function activation is equal to two. 

2.1.2 Tether Bowing Due to the Eleotrodynnmic Drag 

Measurement of the direction of the subsatellito by observation of the 
angle of the tether will be subjeot to error as a result of the curvature 
of the tether caused by drag and eleotrodynamic forces. Calculations have 
been done to estimate the size of the error and provide a possible means of 
correcting measurements for tether curvature. 

A program oalled CURVES is available at SAO for calculating the 
equilibrium shape of a tether under tho influence of eleotrodynamic and 
gravity gradient forces. The program uses the numerical integrator from 
SKYHOOK to generate the shape of the curve given the position, orientation 
and tension of the wire at some point. 

Neglecting the effects of atmospheric drag, the tether must be 
vertical at the point of attachment at the subsatellite. The tension must 
bo equal to the gravity gradient force on the subsatellito. Assuming that 
the ourvature of the tethor is small, the altitude of the subsatellite must 
be that of tho Shuttlo plus the length of the tether. 
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Program CORVES has boon run for the esae of a 550 kg aubaatollite at 
the ond of a 20 ku tether deployed upward from the Shuttle at an altitudo 
of 296 km. The resistance of the tether is .2 ohms/metcr (4000 ohms 
total) , and the current is resistively limited. If the magnetic field (B) 
is .3*10“'* w/m 2 , and the orbital velocity is 7,730 km/sec, the induced 
voltage is about 4600 volts so that the current is 1,15 amps. For an 
orbital angular velocity of .001158 radians per see the gravity gradient 
force on the sobsatellite is 45.8 newtons. The eleotrodynamic force IB is 
.69 newtons. The density of the wire is B.35 kg/km. This information is 
needed to integrate the tension which affects the curvature. The 
integration was started at x ■ 0, snd y » 20,000 meters, where x is 
measured slong trsok and y is along the looal vertical. The generated 
curve intersected the x-azis 133.84 meters from the origin with a tension 
of 52.75 nowtons and an angle of .013080 radians (.7494 degrees). The line 
of sight angle to the subsatellite is .006692 radians (.3834 degrees), so 
that tho error by measuring the tether angle for tracking the subsatellite 
is .3660 degrees (see Figure G-a). At the upper end the deviation from the 
lino of sight is of course .3834 degrees since the integration was started 
with the wire vertical at the subsatellite end. The radius of curvature of 
the wire at the upper end is 1,327,000 ueters and at the Shuttle end is 


1,529,000 meters. 
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Figure G - Geometry for the Tether Bowing 
Analysis (Figure not to Scale) 

(a) - Reference frame and deflection ingles 

as computed by the CURVES compucer code 

(b) - Reference frame and mass discretization 

as in the SKYHOOK computer code. 
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The deviation of the tether angle from tho line of sight con be 
calculated approximately from the curvature of the wire. Tho curvature ia 
related to the tension and the el tfotrodynamic forco per unit length by the 
equation JT7 /d8 = T/R^ , Tho force per unit length ia IB B 3.45 x 10 ^ 
newtons/meter. Tho gravity gradient forco on tho vire is 6.95 newtons and 
that on the subsatellite is 45.8 newtons so that the tension at the Shuttle 
end is 52.75 newtons. The values of 45.8 and 52.75 newtons give curvatures 
of 1,327.000 and 1,529,000 meters as calculated in tho integration program. 
A radius of curvature of 1,529,000 meters gives a change in angle of .75 
degrees ova. 10,000 meters. The deviation from the line of sight would be 
half this or .375 degroes. This is in reasonable agreement with the 
deviation, from line of sight calculated by the integration program. 

A simulation has been run with the SKYHOOK program using 5 mass points 
(three representing tho wire) to further study the deviation from line of 
sight with a current carrying tether (soe Figure G-b) . Tho simulation uses 
a 2.59 mm tether with each mass point having a mass of 41. 7 j kg. The 
Shuttle ma6s is 101.3379 metric tons and the subsatellite is 550 kg. The 
magnetic pole is oriented 90° west of the starting point of the orbit. Ths 
tether is initially vertical and is deflected to the rear, oscillating 
about the equilibrium angle. The program produces a file of in-plane and 
out-of-plane components for each mass point. The pendular part of these 
displacements has been removed in order to study the deflection of the wire 
directly. Figure H and I show the deflection of each mass point vs. time. 
The wire current is also plotted in Fignre J, Mass 1 is the Shuttle, mass 
2 is the subsatellite; mass 3 is next to the subsatellite, mass 4 is in the 
middle of the wire, and mass 5 is next to the Shuttle. The amplitude of 
the in-plane oscillation of masses 3 and 5, soon after the current 
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Figure I - Tether deflection wrt. the Orbiter/Satellite line of sight. 
Pendular libration removed (MKS Units). 
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activation, is about 30 meters and for mas* 4 it is about 44 motors. Since 

each segment of wire is 5 km long, the deflection angle at the ends of the 

wire is about .34 degrees. In order to compare this to the output of 

program CURVES it is neoessary to corrcot for various differences. In 

CURVES the current was 1.15 amps and the magnetic field was .3 x It)** w/m^. 

In the SKYHOOK run the initial current is .81 amps and the z-component of 

the magnetic field is .2 56 x 10 - ^ w/m^. Since the wire is oscillating 

about the equilibrium displacement in the SKYHOOK run, the amplitude should 

be twice as great as the amplitude of tho equilibrium displacement. 

Beoause of discretization the angle computed from the wire mass adjacent to 

the end gives the mean angle about 3/4 of the way from the center of the 

wire rather than the angle at the end of the wire. Assuming the curvature 

is nearly constant, the angle is too small by a factor of about 3/4. 

Taking the defleotion of .38 degrees from program curves and applying all 

the correction factors, 2 x (3/4) x (.81/1.15) x (.256/. 3), gives ,34 in 

agreement with the results from SKYHOOK. In program CURVES the maximum 

deflection at the center (with the pendular part removed) is about 32.5 

meters. To compare this with SKYHOOK tho amplitude must be multiplied by 

the factors 2 x (.81/1.15) x (.256/. 3) - 1.2. Multiplying 32.5 x 1.2 » 39 

meters which compares fairly well with 44 meters from SKYHOOK. The 

amplitude of the oscillations in the out-of-plane direction is about 2 or 3 

"4 2 

meters. The y-componsnt of the magnetic field is .021 x 10 w/m which is 
8% of the z-componcnt. Multiplying 30 meters by .08 gives 2.45 meters 
which agrees with the observed amplitude of the out-of-plane oscillation. 

In addition to tho oscillations, the in-plane and out-of-plane displacement 

-V 4 

show shifts in the mean value through the orbit as the magnitude of the IxB 
force slowly changes. The period of the transverse oscillations is about 
535 seconds both in-plane and out-of-plane. The period of the fundamental 
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fundamental transverse mode Is 2 2/v where the velocity of propagation v is 
^T/p . Taking T = 46 nowtons, p ** 8.35 x 10~^ kg/m, and £ D 20 km gives a 
period of 539 soconds in agreement with tho observed period. 


2.1.3 Concluding Remarks 

The satellite and tether dynamios during the entire station-keeping 
phase of the first eleotrodynamic mission does not appear to be worrisome. 
The out-of-plane satellite libration builds up to a maximum value of 1.3 
deg in a worst condition case: current switched on abruptly, plasma 

contactors at both tether ends, tether resistance of 4000 ohm (as presently 
planned by MMA) , and no other resistive load in the loop. The build-up is 
due to a small resonant component (frequency = 20) of the eleotrodynamic 
force. The system is however forced out-of-resonanoe by (most probably) 
the magnetic pole rotation. Damping out tho out-of-plane swing amplitude 
by tether tension control came out to be uneffective. Current control (if 
deemed necessary) could be a valid option. 

The tether bowing is fairly limited. A conservative value is obtained 
by multiplying by the dynamic overshoot (a factor of 2) the result obtained 
with the program CORVES. The maximum deflection (in the middle of the 
tether) from the line of sight is, therefore, around 65 m. This value is 
consistent with a maximum angle (at the boom tip) between the tether and 
the line of sight of 0.767 degrees. 
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2.2 High Resolution Slack Tether Studies: SLACK2 Model* 

2.2,1 Introductory Remarks 

Under previous contract, SAO has developed a high-resolution, two 
dimensional model of a zero-tension tether (Gullahorn, 1983; Colombo, 
Arnold, Gullahorn and Taylor, 1984). 

This model, implemented in the computer program SLACK2, is intended to 
be similar to SKYHOOK insofar as feasible. As with SKYHOOK, the continuous 
tether is modelled with tho simple "lumped mass" or "ball and spring" 
approaoh. A new physical system is constructed by dividing the tether into 
several segments; each segment is replaced by an equivalent mass at one 
end and a massless spring connecting this mass to the adjaoont segment's 
mass. The spring generates no force when not stretched; the segment area 
(for air drag) is mapped onto the end mass. SKYHOOK then applies realistio 
forces, both internal and external, to each mass and numerically integrates 
their motions in earth-centered coordinates. SLACK2 operates in a 
coordinate frame relative to the (infinite mass) Shuttle and applies an 
idealized set of forces: linearized gravity gradient and Coriolis 

accelerations; drag; perfectly inextensible (infinite elasticity) 
connecting segments. This latter idealization means that the segments are 
"almost always" slack and that when two masses separate to bring a segment 
into tension they undergo an infinitesimally short "bounce" reversing their 
motions along the segment direction in their center of mass frame. Both 
the bounces and the motion of the masses between bounces are analytically 
soluble, allowing a major increase in efficiency over SKYHOOK, which must 
numerically integrate difficult differential equations. The motion is 
restricted to the orbital plane in SLACK2, but this is not a fundamental 
requirement and a three dimensional version is in preparation. 

★ 

Contributed by Dr. Gordon E. Gullahorn, SAO 
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la contrast to these restrictions for improved speed, SLACK2 eztonds 
SKYHOOK'S capabilities in other directions: initial conditions are 

generated internally, and non-uniform segment lengths are supported; many 
more segments aro feasible, up to fifty; a vibrating boom is included as 
the tethor attachment point. For more details on the model, seo Colombo ot 
al. (1984). 

In tho current reporting period, SAO has improved and studied SLAC&2 
and utilized the program in tether break studies: 

— The post-processing plotter was re-created for operation with 
a changed computor operation system. A printer plotter/lister 
was oreatod to allow detailed examination of particular variables. 

— Improvements in the coding and some of the numerical routines 
were made which reduce the run time by a factor of two. 

— A detailed comparison was made between SKYHOOK and SLACK2 results 
for a particular case. All discrepancies are understood. 

— A protoool for running cases with long tether remnants was 
developed. Previously, such cases had caused program failure. 

— The effects of tether longitudinal damping were studied. 

— Modifications toward a threo dimensional version are underway; 
initial coding of a simple version is complete. 

— A stability analysis of a very simple version (linoar, 
no Coriolis or drag forces) was initiated. 

Several of these topics are discussed in further dotail in the 
sections below. 

2.2.2 Comparison of SKYHOOK and SLACK2 results. 

SKYHOOK has been exercised, debugged and refined for many years. 

Since the physical model in SLACK2 is similar to that in SKYHOOK, an 
. obvious way to enhance confidence in SLACK2 is to oompare its results with 

those of SKYHOOK. For a variety of reasons, making this comparison is not 
as trivial as one might at first expect: 

w- — ■ ■ . .. • 
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— SLACK2 incorporates an oscillating twenty meter boom; 
this is not in SKYHOOK. 

— The initial conditions for SKYHOOK are tedious to oompntc (by hand) 
unless one wants to start from an equilibrium configuration with 
equally spaced masses, 

— SKYHOOK assumes a tether with finite elasticity; SLACK2, infinite. 
Thus, the instantaneous "bounoes" in SLACK2 become brief time 
periods in SKYHOOK during which tho motion is integrated with a 
stretched tether. Theso intervals become briefer, and tho agreement 
with SLACK2 better, as the elasticity in SKYHOOK is increased, but 
thero is a limit beyond which SKYHOOK experiences difficulty. 

— SLACK2 assumes a constant density, and this constant is built into 
the program (for ease of operation). SKYHOOK oomputes tho density 
as a function of altitude and solar direction, 

— SKYHOOK computes the mass and area/mnss, rather than accepting these 
as input, for eaoh tether segment ... except the terminating 

one. This exception is undocumented and not clear from the code. 

The values for the terminating segment must be input as the 
"subsatellite." (If we were operating with an actual subsatellite 
instead of a cut tether, the subsatellite mass and area/mass would 
have to ohange with the number of tether segments to be correct.) 

— SKYHOOK and SLACK2 use somewhat different drag models: in SLACK2 

area of each segment never varies; in SKYHOOK, the area 
perpendicular to the direction of motion is used. SLACK2 uses 
the total segment area, multiplied by a fudge factor of 0.75 to 
approximate the projection effect. The SLACK2 model was chosen 
primarily for computational reasons, but it is not a priori 
obvious that the SKYHOOK projected area model is more appropriate 
for a slack, possibly crumpled, tether segment. 

— SKYHOOK and SLACK2 use different conventions for the sign of the 

"in-plane” coordinate in tho plot file: in SLACK2 the "in-plane" 

axis points in the direction of orbital motion (a positive value 
indicating further "along orbit" than tho reference Shuttle), 
while SKYHOOK chooses tho axis pointing back along the orbit 
(this was found by experiment, and was apparently done to make 
the successive, configuration plots more consistent since the 
Shuttle moves from left to right in the successive plots and 
with tho normal conventions the positive in-plane direction 
would be to the left). 


Keeping these considerations in mind, the comparison runs were made 
(after several iterations) as follows: 

Step 1*. SLACK2 was run using the VAX Debug facility. This allowed us 


to halt the program and 




r-v 


^ -.ft SSR ?■».-**> -:••*: - •*»**»*•, V. twa^ V W2f T7 ±K^AW'ES^5* l V 1 


Page 31 


— act the boom length to zero, eliminating the boom effects, 

— set the atmospheric density to tho value found from previous 
SKYHOOK runs at the given altitude, 1,91 x 10 1 gm/cm , 

— set the drng "fudge factor" to 1.0 
beforo the aotual simulation. 


Step 2: An input file for SKYHOOK was created and subsequently run. 

This was done by 


— starting with a typical equilibrium configuration input file 
as a model, 

— obtaining the initial vertical displacements and velocities 

of the "masses” relative to the Shuttle from the SLACK2 output 
and from internal SLACK2 variables examined with the debugger, 

— converting these to SKYHOOK initial conditions by 

- adding the displacements to the Shuttle vertical coordinate 

- soaling the horizontal velocities to give the same angular 
velocity as the Shuttle 

- putting in the vertical velocity 

— setting the Shuttle mass to a high value (10^ gm) and 

area/mass low (0.36 x 10 °) to approximate idealization 

in SLACK2 , 

— setting the "subsatellite" mass and area/ma66 to the values for 
the final tether segment, 

• — putting in a high tether elasticity, E 0.7 x lO*^, about 

a factor of 100 higher than Kevlar, to approximate the infinite 
E idealization of SLACK2 (higher values of E lead to 
SKYHOOK failure). 


The case run in both programs corresponded to the following parameters 
(they were selected for this comparison only): 


— altitude 220 bn, circular orbit 

— 0.2 cm diameter Kevlar tether 

— original 100 fan tether deployed upwards with 0.3 ton 
subsatellite (determines recoil velooity) 

— tether broken 0.2 fan from deployer 

— discretized to five tether segments (6 masses in SKYHOOK, 
including the Shuttle; 7 in SLACK2, with Shuttle and boom tip) 

— evolution followed for about BO seconds. 
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Comparing and analyzing the resulting outputs is complicated by a 
variety of factors: the above mentioned disagreement in sign of the 
in-plano component; the presence of tho deployment boom tip as an extra 
"mass" in the 3LACK2 output; lack of a standard metric for comparing 
tether configurations. Development of the lattor (perhaps a sum of squares 
of differences in mass position and/or differences in angles betwoen 
segments) would be valuable for stability studies; but more intuitive 
considerations were adequate for the current comparisons, and highlighted 
many of tho problems which had to be dealt with in early iterations. Three 
factors were considered separately: the gonoral appearance of the 

configuration (side-view) plots; the vertical positions as a function of 
time; and tho horizontal in-plane accelerations. 

In-nlane accelerations : This factor proved to be the most immediately 

diagnostic of differences between tho models, leading to the 
corrections/allowances detailed above. The in-plane positions of the two 
end masses were extracted from the plot file using the printer 
plottor/lister. The resulting list of position as a function of timo was 
then differenced twice to give acceleration (and the negative taken of 
SKYHOOK results). 

In the initial period before any interactions, we expect a constant 
acceleration due to Coriolis forces and differential drag. This was indeed 
observed, and ufter discrepancies were identified and compensated for, the 
two programs gave identical accelerations: 1.20 cm/sec^ in tho forward 


direotion 
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Time history si vortical coordinates i Tho vertioal ("radial") 
coordinates of the mosses with respect to tho Shuttle need no adjustment of 
sign and were plotted with the standard SKYHOOK poet-prooeesor for the 
results of both programs. Those are shown in Figure K-l, Note that the 
SLACK2 results includo a seventh masB, tho boom tip, which is coinoidont 
with the Shuttle and hence hus constant coordinate 0. The results are 
almost identical up until about T - 40 sect after that they bogin to 
diverge noticeably. Doteiled examination of tho printed results shows that 
this divergence coincides with a pair of nearly simultaneous bounces 
involving masses 3 and 2, then masses 2 and 1, between T *» 40 and T 13 41. 
Tho bounoe times are printed explicitly by SLACK2 but mu6t be inferred from 
the history of integration steps in tho SKYHOOK output, the integrator 
taking very small steps near any discontinuity such as a change from slack 
to tensioned for any segment. For lower values of elasticity in SKYHOOK, 
these two bouncos were reversed in order and the subsequent divergence was 
much more pronounced. 


Configuration nlots : A complete set of successive configurations at 

one second intervals is shown in the familiar in-plsne vs. radial 
component side-view plot for the SKYHOOK results in Figure K-2. To allow a 
more detailed comparison of rosults we have plotted side-by-side the 
configurations of the SKYHOOK and SLACK2 results for several values of T in 
Figure K-3, oppropriately inverting the SKYHOOK configurations to 
compensate for the different in-plane sign conventions. Results are in 
detailed agreement up to T «= 40; reasonably good at T = 50; and agree 
only in a very general way thereafter. Of course, these are comparisons at 
the same t jae and do not rule out tho possibility that the primary 
difference is simply a change or jump in tho time scale, as can happen with 
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orbits: two particles in very olose orbits will drift far apart if looked 

at simultaneously. 

In summary, we may say that the agreement between SLACK2 and SKYHOOK 
is excellent once known differences are compensated for. Thoso differences 
were intentionally chosen, in some oases representing extensions or 
corrections to SKYHOOK, in some cases for efficiency. The divergence in 
results after T « 40 sec is due more to the limitations of the overall 
"Ball and Spring" model than to problems with the SLACK2 program: tho 

different results of the series of bounces between T = 40 and T = 41 could 
as easily have resulted from very slightly different initial conditions as 
from the slightly different treatment of bonnoes. This raises questions 
about the stability of the "Ball and Spring" modol, and indeed of the 
actual physioal system, which have only begun to be addressed. 


2.2.3 Long Tether Remnants 

In the SLACK2 model ell free-flight forces are linear. Thus, when we 
start from the basic configuration of a set of masses traveling at the same 
velocity and separated by tether segments each slackened to the same degree 
(i.e. with separation the same fraction of natural length), all segments 
will come into tension simultaneously. When only a short length of tether 
remains after the break, interaction with the boom (to which the remnant is 
still affixed) disturbs the tether before the uniform lengthening can cause 
problems. Longer remnants, however, become taut before even the first mass 
is affected by the boom. Although SLACK2 can handle occasional nearly 
coincident bounces, there is no obvious algorithm for handling a case where 
several occur so as to interact: the program usually enters a broadly 

defined loop, bouncing the masses in some sequence over and over, or simply 


fails with some fatal error. 


The method used to overcome this problem is to make the initial series 
of bounces non-simul taneous. This is done by introducing some 
randomization in the initial conditions. It is found experimentally that a 
few percent randomization in the slackness factor [separation/ (natural 
length)] used for different segments, and one or two degrees randomization 
in the direction of the segments, is sufficient to prevent program failure. 
Figure K-4 shows results for a 100 km tether cut at 20 km, a case which 
previously could not be handled. 

Typically there is a sharp buret of bounoes when the tether initially 
tries to come into tension, the masses sort out their motions, and settle 
down to a regime of more normal bounce rates. E.g., in a case similar to 
that illustrated, there were no boun'es until T =66 sec, some 25 bounces 
until T = 69, 369 between T =69 and 70, and thereafter a slow oscillation 
between a few bounces per second and a few tens of bounces per second. 

2.2.4 Stability of the Lumped Mass Model 

This resolution of our operational problem introduces seriously the 
questions about stability. In order to run at all, we must perturb the 
problem slightly; but do all such perturbations yield some identifiably 
similar end result, at least in the immediate aftermath of the initial 
burst of bounces? In general, for how long are simulations beginning with 
slightly perturbed results only slightly different? Those questions also 
arose in the comparison of SKYHOOK and SLACK2 results, where the 
perturbation was to the details of the model rather than to initial 


conditions. 
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Results from a SLACK2 run with a long tether 
remnant after break (20 kilometers remaining 
from 100 km tether; 30 segments in simulation) 
This demonstrates that a small randomization 
of the initial conditions allows successful 
simulation in long remnant cases. Configura- 
tions are shorn at 50 second intervals for 
2500 seconds. 
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Investigations of a one-dimensional model, with no attachment point 
and only the gravity gradient force, were begun. Even a simple three-mass 
model leads to rapidly expanding numbers of special cases when studied 
analytically, before any significant results are obtained. We hope to 
program this modol, perhaps on a microcomputer, and oxamine results of 
varying initial conditions. 

2.2.5 Tether Longitudinal Damping 

In this reporting period we used a crude damping feature that had been 
built into SLACK2 but not previously exercised. 

Consider two masses and the connecting spring in their center of mass 
frame. If we assume a typical damping in the spring proportional to the 
rate of stretching, then each "bounce" {the one-sided spring becoming taut, 
stretching from the inertia of the masses, rebounding and becoming slack 
again) will result in a decrease of the masses' relative velocity by a 
constant fraction. This featuro has been programmed into SLACK2 in that 
one inputs a "percent damping", e.g. 2% results in post— bounce velooity 
0.98 times the pre-bounce velocity. This damping factor is applied 
uniformly to all segments regardless of length. (In the future we hope to 
investigate the appropriate physics so that we can input a physical 
quantity related to the tether material, and compute the appropriate 
damping factor for each segment.) 

We shall use below a quantity we call the "retention". This is a 
simple approximation to the velooity we would expect to remain after a run, 
and i6 computed (post-hoc) as follows: count the number of bounces and 

divide by the number of segments to get the number of bounces per segment; 
if a segment were bouncing in isolation, after N bounces its velocity would 
be (damping factor)^ times the original, so we define this retention as 
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(damping faotor) I( bounoC8 > /( segments)] 

A sequence of runs with different damping factors were made for the 

case : 


— original 20 km tether deployed upward with 
0.5 ton subsatellite 

— cut at 0.2 km 

— boom deployed forward 30 degrees 

— discretized by 35 segments 

— evolution followed for 600 seconds. 

The results are shown in Figure K~5 for damping factors of 0%, 0.5%, 1% and 
2%; the resulting retentions were 1.00, 0.45, 0.14 and 0,02. The results 
of increased damping are about what one might expect: there is an overall 

decrease in the "crinkly" nature of the configurations and the effect of 
the vibrating boom becomes less apparent; the bounoing of the tether as a 
whole about the attachment point has longer period; and the effects of 
drag become more pronounoed. An unexpected feature is that the first 
excursion past the attachment point is further the greater the damping; 
the boom seems to retain some tether in its vicinity in the undamped case. 

The 2% damping case was followed for a longer period, 3000 seconds. 

The results are shown in Figure K-6. Note that the tether appears to be 
oscillating nbout a "downwind" dragging position. With the "balloon" drag 
model used, and approximating the tether as a single mass, this downwind 
configuration can be shown to be stable for tethers less than about a 

kilometer long. The retention factor for this run is extremely small: 
(0.98) C<19058 bounoes)/(35 segments)] = % % 1Q -5 # 0np wonld expect t0 8ec 

a perfectly motionless system, but instead the tether is obviously still 
moving moderately briskly at least. Possible sources of energy input are 



f 
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air drag and the vibrating boom. It would be of interest to compote the 
total kinetio energy as a function of time; we might e.tpeot this to show 
an initial decrease duo to damping but eventually reach an equilibrium 
between damping and forcing. 

We should point out that under the influence of damping the ''almost 
always slack" model used must eventually break down. In a real tether, the 
gravity gradient and/or drag forces will eventually bring the tether into 
constant tension, unless there is some disturbing influence such as tho 
vibrating boom. How will SLACK2 respond to this challenge? A run with a 
non-vibrating boom should be made, but one may make a thooietioul 
prediction as well: If one analyzes the simplest case cf one tether mass 

and an infinite mass Shuttle with only the gravity gradient force applying, 
one finds that the velocity reduction at each bounce causes the time to the 
next bounce to be reduced by the same ratio. Thus, as the system 
approaches tension, tho bounces become infinitesimally spaced and approach 
a constant time. We have not seen this directly in the SLACK2 runs because 
of the booms energy input. 


2.2.6 References 


Colombo, G. , Arnold, D.A. , Gullahorn, G.E. and Taylor, R. S. , 1984. 
"Investigation of Bleotrodynamic Stabilization and Control 
of Long Orbiting Tethers," Final Report on Contract NAS8-35036. 

Gullahorn, G.E., 1983. "Slack Tether: High Resolution, Two 

Dimensional Model," technical appendix to Monthly Progress Report 
#11, Contract NAS8-35036 (G. Colombo, PI). 
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2.4, Preliminary Analysis of Eleotro-Ueobanloal Failares 
Cid Tentative Asaeaament of the Related Consequenoes 

2.3.1 Introduction 


Tho elect rodynam ic tethers that will be flown in the forthcoming first 
and third demonstration flights of the Shuttle-borne TSS facility oould, in 
principle, represent safety hazards due to tha high-voltage electro-motive 
force that they will gonerate. In this preliminary analysis we show that 
these hazards ore not too worrisome. Thus far. however, we have limited 
our investigation to tho hazards to b6 expected from possible failures 
occurring while tho 20 km tether i6 fully deployed (upwards). Bofore tho 
oocurronce of failures, the lower end of the tether is assumed to be 
connected to the Shuttle and the Shuttle is assumed to bo kept at tho 
plasma potential by a suitable plasma contactor (hollow-cathode plasma 
bridge). The upper end of the tether is assumed connected to the 
subsatellite. Two different situations are assumed for the latter: (1) 

presence, and (2) absence of a plasma contactor. 

The main failure modes are the following: 

1) Malfunction of the Shuttle-borne plasma contactors; 

2) Malfunction of the plasma contactor in tho subsatellite (when such 
a plasma device is present); 

3) Breakage of the tether (either near the upper end, or in the 
middle, or near the lower end). 


The simplified diagram of the eleotrodynamic tether system is given in 
Figure L. For a 20 km tether, the electromotive force (emf) hns an orbital 
average of 3.78 KV, for a typical 28° orbital inclination and a 295 km 
orbital height. The contact resistances R^ and are of the order of IQ, 
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wbon plasma oontaotors arc mod at both tether's ends, and when they 
function proporly, The roelstanoo of tho ionospheric olosuro circuit Rj ono 
is also of tho ordor of 1 ohm. Prosont plans for the 20 km wire oall 

cither for R^ » 1,686 kfl (25.7 ohm/1000 ft) or for 4k Q (.2 ohm/m). Wo 
assume here that is also equal to either 1,686 Q, or 4000 Q. 



Jf.M t '.t 
( -t Air 

Fig L Simplified equivalent citrun of elcctrodj namic tether (tether deployed 
upwards, shuttle in eastward motiunt. 

Consequently, for a 1.686 kohm tether the current 1^ in the wire is 
expected to be about = 1,12 A, with 2118,6 watts di««ipatcd in the 
tethers ohmic resistance and an equal 2118.6 watts dissipated in the 
on-board loads. In fact, when plasma contactors are used at both ends of 
the tether, a current intensity « 1.12 A is well with feasibility. On 
the contrary, when plasma contactors are used solely on the Shuttle 
Orbiter, tho current 1^ in tho wire is limited by the charge-capture 
capability of the terminating subsatellite (a metal sphere with diameter <= 
1.5 m) . This upper limit for the current has been computed to be ■= 0.63 


A. 
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2.3 .2 Consequences of Plasma Contactor Failure 

If plasma contactors ore used at both ends of the tether, failure of 
one of them will have a 3 a consequence the appearance of a high potential 
between that termination where the hollow-cathode plasma contactor failed, 
and the surrounding medium. Let's consider the case of hollow-cathode 
failure on-board the tethored subsatellite. When thin occurs, the 
"contact" with the ionosphere is provided now solely by the metol skin of 
the 1.5 a diameter spherical subsa tel 1 ite. We have already indicated that 
the current will drop from 1.12 A to 0.63 A. The resistance of the 
equivalent circuit (see Figure L) has now increased in value from about 1 
ohm to 2,628 Q. A difference of potential + 0.63 * 2628 «* +1665.6 volt is 
established between the external surface of the eubsatellite and the plasma 
environment. This potential is highly super thermal and ve should not 
expect that its distribution around the subsatellite will follow the 
prediction of quiescent plasma sheath theory. It is more likely that, 
owing to collective effects, the perturbed volume will be larger than 
predicted and will resemble some cort of "corona discharge" with 
accelerated plasma populations inside. Hazards, deriving from this 
situation do not seem, however, to be worrisome. It is more probable that 
we will have here the opportunity of performing interesting plasma physics 
measurements, rather than experiencing damages to the subsatellite 
instrumentation. The risk at least seem affordable. Instead, a serious 
hazard is represented by the possibility thnt the Shuttle-based 
hollow-cathodes may fail, thus leaving the electric potential of the 
Shuttle free to inorease from such a safe operational values as a few volts 
(negative) to very high negative values with respuct to the surrounding 
medium. An initial calculation (to be verified in next reporting period) 


V ^wKCSt-W' *' * 1 *!-■*»: :,„***«► , #*a.in * **■ 


Page 50 


performed by assuming that the upper end of the tethor is terminated by a 
1.5ra diameter subsatellite without plasma contactor, has indicated that the 
negative voltage acquired by the Shuttle in case of failure of the 
Shuttle-mounted hollow-cathode is about -2kv with respect to the 
environment plasmn. 

We plan to perform in the forthcoming contract's performance period 
detailed calculations of the potentials acquired by each tether's end in 
cose of failure of plasma coutactors. It is already obvious, however, that 
a hollow-cathode failure at the Shuttle end would have serious 
consequence s. Provisions must be taken (by use of a redundant arrangement 
of several hollow-cathode devices in parallel, or by some other approach) 
to reduce as much as possible the probability of the occurrence of such an 
event. It would also appear that, in order to alleviate at least in part 
the consequence of the failure of the Shuttle-mounted hollow-cathodes, the 
plasma contactor on-board the subsatellite should be switched off by the 
operator. This would make it possible to localize part of the excess 
difference of potential at the subsatellite end of the tether (where it is 
less worrisome) thus reducing the negative voltage acquired by the Shuttle 
Orbiter with respect to the plasma medium, 
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2.3.3 Consequences of Wire Breakage 
2. 3. 3.1 General 

In all the analytical efforts that have .. . u performed thus far on the 
eleotrodynamic interactions between the ionospheric medium and an insulated 
metallic tether with terminating electrodes, the assumption always has been 
that the typical dimension of such electrodes are greater ihan the electron 
Debye length At the orbital altitude of the TSS demonstration 

flights, ■= 3 to 5 mm. This condition is amply met when the free end of 
the wire is terminated by a subsa tel 1 ite . This is however no more so, when <| 

a wire breakage occurs. In this case a high potential could appear betweon i 

the small dimension, truncated tether's tip and the surrounding medium. A 
first-out estimate of the extent of this high potential was performed by 

Olbert (1983), who modellod the long tether as an elongated prolate j 

I 

spheroid in vacuo. Olbert' s approach would lead to values of the eleotric j 

♦ 

field near to the truncated tip of the wire that would be high enough to be j 

i 

i 

worrisome. However Harrold (1984) has shown that we can realistically \ 

expect there lower values for tho electric field. 

, \ 

Section 2.J.3. 2 of this report illustrates Olbert's approach, while 
Section 2. 3. 3. 3 contains Harrold's revision and his estimate of the i 

expected values of the eleotric field near the tip of the truncated wire, j 

i' 

at various distances from the wire's tip. 1 
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2.3 .3.2 Olbert's Method of Computation and Numerical Results 

As already pointed out in the previous Section, Olbert (1983) modelled 
the long wire as an elongated prolate spheroid vacuo. Qis computation of 
the electric field noar the truncated tip of the wire proceeds as follows. 
With reference to Figure M, the electric field E^ at the point of the wire 
that is at a distance y from the tether's mid point is given by the 
expression 

y a 

00 

I E | — (1 ) 

Qjb/a“- (a 2 -b 2 )y 2 


where (see Figure M for definition of a, b): 

<= field at infinity ° |V x S| 
a = tether length (after breakage) 
b » wire's radius <= 10 J m 

°x - >" < t ! -*• 


The highest value E occurs when y = a. in this case, equation (1) 


yields : 


(max) 


(max) 


2 „ 
a E 

00 

~2~~ 
b Qj 


( 2 ) 


Because, in our case E m = V x S| •= 0.189 V/m, wo obtain from (2) 


n 


(max) b Q 


a V 

^ x 0. 189 - 

^ - m 
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Infinite 

Ground plane 
(at the Shuttle) 



a = semimajor axis of equivalent 
prolate spheroid = ^xlO^m 


b = Semiminor axis 


In-vacuo, elongated prolate spheroid 
model of long orbiting tether 


' SFV^vazf^’**-** ***•«*• •**•**•<■ 

WE;, k ^.’ , .' i 


Page 54 


Lot's examine the following oases: 


1. Wire breakage at attachment point to subsatellite (20 km above 
Shuttle) : 

a ■= 2 x 10 4 m 
<*! “ 16.5 


‘(max) 


4*10 


E_ = 4.5 x 10*^ V/m (at truncated end) 


10“ 6 x 16.5 


2, Wire breakage, 2 km above the Shuttle: 


a “ 2 x 10^ m 


°1 = 14.2 


E 


“(max) = E ” “ 5.3 x 10 10 V/m (at truncated end) 


10 u x 14.2 
3. Wire breakage, 200 m above the Shuttle: 
a - 2 x 10 2 m 
^ « 11.899 


4x10 


,8 


”(max) 10" 6 x 11.899 
4. Wire breakage, 20 m above Shuttle: 
a = 20m 
^ ■ 9.5966 


E^ ** 6.3 x 10 V/m (at truncated end) 


E 


4x10* 


“(max) 10" 6 x 9.5966 


E„ «= 7.9 * 10^ V/m (at truncated end) 


5. Wire breakage, 2 m above Shuttle: 
a = 2m 


= 7.29 


“(max) 10 6 x 7.29 


E^ *= 1.03 x 10 3 V/m (at truncated end) 



#** 
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2. 3. 3. 3 Revision of Olbert's Approach* 

2. 3. 3. 3.1 Introductory Remarks 

The eleotric field in the vioinity of a wire travelling through a 
uniform magnetic field B can be calculated from the approximation that at 
velocities much smaller than the velocity of light* the Lorentz 
transformation states that the wire is immersed in a uniform eleotric field 

fi^xg (3) 

As viewed by an observer travelling with the wire. 

One approach to the calculation of the distortion of tho otherwise 
uniform E field is to model the wire as a long slender prolate spheroid 
(Olbert, 1983). It will be shown here that while this modol is reasonable 
at some distance from the ends of the wire it is very inaccurate close to 
those ends. 

The electric field calculated at the end of the wire, by the Olbert's 
method is on the order of a teravolt/meter. Even if this erroneous result 
were correct, it will be shown in this section that the electric field, one 
meter from the end of the wire, would already be down to 114.64 
volts/meter. 

The assumptions made for the calculations carried out in this section 
are as follows: 


•Contributed by W.J. Harrold (see reference Harrold, 1984, in Section 

2 , 3 . 3 . 3 . 6 ) 
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a - length of wire 20 tan 
b - diameter of wire 2 mm 

o - one end of wire attached to a largo conducting 
spacecraft (Shuttle Orbiter) 
d - prolate spheroid model 

— semi major-axis a «* 2 x 10^ m 
—semi minor-axis b « 1 mm 
e - spacecraft modelled as an infinite 

ground plane so that the tether length 
is doubled and therefore a = 2 x 10^ m 


2* 3. 3. 3. 2 Prolate Spheroid Model 

It can be shown that the electric field on the axis of a dielectrio 
spheroid, of dielectrio constant K, immersed in a uniform field E o is given 
by 


where : 


C 1 [Q o ( ^ " -e— 1 “ E 


c 2 -i 


■'o' 


z>a 


(4) 


? = z/o 


(5) 






*■ *nf«^ **8S> ***e**1PW£wBP tr - *»*»■*• ' tw**** 
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z *• the distanoe from tho center of the spheroid (meters) 

2 2 2 

C •* the semi-fooal length of the spheroid C ■ a - b 
a “ the semi-major axis 
b ■ the semi-minor axis 


i. in ( J^ti) 

2 £-1 


( 6 ) 


and : 


C ? oEo 

1 = ^ (C c ) + l/tfk-lX^-D] 


(7) 


where : 


K 0 = a/c 


( 8 ) 


and : 


Q 1 (C) = « 0 (U-1 (9) 

The function Q o (£) and Q o (£) are the zeroth and first degree legendre 
functions of the second hind of order o. 

One can calculate the field around a conducting spheroid by taking the 
limit is k + “ thus: 


Lim C 

!(-*■ » 1 


€ 0 Eq 

Q l^o ) 


( 10 ) 


■5 *o£* 




2. 3. 3. 3. 3 Eleotrio Field at the End of the Spheroid 


We substitute (10) into (4) for £ » K and find that 

o 

E n = oflf^T 1 ~ E 0 


E 


[<MO + 1 


Q^) LVi l^o 
E 


... ] ’ E o 




QiC« 0 ) " C 2_! ] " E c 


-E. 


(es-UQiCCo) 


An approximate calculation of can be performed as follows: 


V 1 


where: 


Therefore: 


W = f^ (fry) - 1 


c 0 e a /0 
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( 11 ) 


( 12 ) 


(13) 







W - f ln <srf > - 1 


(14) 


where: 


= - b 2 

- w 


(15) 

(16) 


~ = 10" 3 / 2 x 10 4 = 5 x 10' 8 


and: 
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Therefore C is closely approximated by: 


o « * [1 - 


I <;> 2 ’ 


a - 


hi 

2a 


and 


a - c 


hi 

2a 


Thus : 


i /• a+c '> 

In ( ) 

v a~c 


■ In ( &k } 

’ 21 " <tt) 


Substituting this result into (13) we find that; 


«1 «o> - j 1- <£> - 1 » 1 . <$> - 1 


Furthermore wo find that: 


e 2 „ - i - h - i » ^ 


We now substitute (19) and (20) into (11) with the result: 

a 2 I E„ I 


E 


n. 


b 2 [ln(^) - 1] 


4.5 x 10 12 V/m 


in agreement with Olbert (1983) [see equation (2) in Section 3.2] 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 
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n 


i 

2, 3. 3. 3. 4 Calculation of the Eleotrio Field on the Axis 
of the Prolate Sphoroid 


The results of tho previous section would like worrisome, if wc do not 
check how small is the region where tho electric field is so high. A 
oomputor program was writton to calculate as a function of (x - t), A 
summary of these results is given in Table I, herennder (see Figure N) . 


Table I 


Values of the Electric Fiold on the axis of tho prolate 
Spheroid at Various Distances from the Tip 


(z - a) 

1 

E ! V/m 
z 1 

10~ 10 m (1 A) 

1.1451 

x 10 12 

l(f 9 ti 

(10 A) 

1.1451 

x lO 11 

10“ 8 m 

(100 A) 

1.1451 

x 10 10 

10“ 7 m 

(1000 A) 

1.1451 

x 10 9 

10“ 6 m 

(l|i) 

1.1451 

x 10 8 

10~ 5 m 

(10(1) 

1.1451 

x 10 7 

10" 4 m 

(100(1) 

1.1452 

x 10 6 

10“ 3 m 

(1 mm) 

1.1452 

x 10 5 

10“ 2 m 

(1 cm) 

1.1453 

x 10 4 

10 _1 m 

(10 cm) 

1.1453 

x 10 3 

1 

m 

114.64 

10 

m 


11.59 

100 

m 


1.30 

1 

km 


0.28 






m 
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2. 3. 3. 3. 5 Demonstration that the Near-Pield of tho Prolate 
Spheroid is an Inoorroot Model 


Tho radius of curvuturo of a curve la given by the reciprocal of the 
second derivative. For example in the oasc of the circle: 


2 2 
x + y 


“ R 


2 


y = 




n . ,-: X .. . 

dx 



dx 2 


-1 

/R 2 -x 2 


CR 2 -x 2 ) 3/? ’ 


At x « 0 : 


a~2 R 


The spheroid nnalyzed in this section is a result of rotating tho ellipse 
defined by: 




1 


About its major axis "2a" 


x 



y 2 /b 2 


dx 
d y 




-a 

b 2 /l - y 2 /b 2 


+ 


(1 

b 4 



- 3/2 


d 2 x 

dy 2 




r * ■■-W'i W>irVW TWJ T»’ 
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At y “ 0: 


or : 


£l . _ _2_ 

dy' ? ’ b 2 



5 x lO" 11 a 


Since the radios of the wire is b we find that the radius of the sharp end 

»"8 

is too small by a factor of b/a or 5 x 10 

The voltage at the surface of a charged sphere can be computed as 
follows. The electric field of a charge q is given by 


E 


r 


q 


4irr 2 Ey 


The potential at a radius r, is given by: 


Thus : 




“ c[ dr 
4rrr 2 E 0 


a 

4 «iE 0 


1 " 4nr l V rl E 0 


and: 



Vrl 

r l 


n 

If r^ is erroneously calculated to be too small by the factor 2 x 10 then 
the electric field at that miniscule radius will be too high by the factor 


3 ■' ■ ^ J 7~» ^ *"• '* ■ . rn.- ' w-» • ■ -■-“ '”■■*'■* *?•" -w.t" : - 5.^ v v --JK » J 7* MV. &£. f Mill- uU *?.n s* 1 

V* * V 

4, * *t*A 
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2. 3. 3. 3. 6 Concluding Remarks 

We have calculated the electric field at the end of a prolate spheroid 
moving in a uniform magnetic field and we have shown that this model 
{Olbert, 1983) gives electric fields which are too high by a factor of 20 
million. We have also shown that, even if one aocepts the predictions of 
the spheroidal model, the electric field strength drops to 114.64 
volts/meter at a distance of one meter from the end of tho spheroid. 

Thus far, our analysis has not disclosed eleotric hazards that are 
serious enough to require some ohangos in NASA plans for the two 
eleotrodynamic tether missions of 1987-1990. We will continue, though, in 
our search, especially in connection with slaok-tether situations that 
might bring the truncated tip of the wire (with corona discharges around 
it) close to the Shuttle Orbiter. 
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2.4 Program Status of the TSS Rotational Dynamics Hodcl 

The modelling of the external perturbations affecting the rotational 
dynamics of the satellite has boon started. The implementation in the 
computer code has not yet been done so that presently tho crude original 
version is the only one running. With regards to the external 
perturbations this version models the tension only, and it is therefore 
suitable for first approximation simulation. We are also working on tb? 
analytical formulation of the rigid body motion that presently integrates 
the nine directions cosine along with the three rigid body equations of the 
satellite. The idea is to derive the expressions of the satellite 
rotational velocities, P, Q, R (in body axis) as a direct function of the 
Euler's rotational velocities. This strategy will decrease the total 
number of equations for tho rigid body rotation of the satellite from 12 to 
6, The disadvantage is to have a singularity at tho equation describing 
tho roll motion for a roll angle of 90 degree from the local vertical. 

Such value however represents a limit condition. First results of the 
up-dated rotational dynamics model will be presented in the next quarterly 
report. 


2.5 Problems Encountered During Reporting Period 


None 




i. y.xynti'.'-: '■=■ ti « ; »v.r,j.ir ; i r," a $*? ./*«iy? 'R»w r ^ #&***?? •*««** 
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2.6 Activity Planned for the Next Reporting Period 

The activity on tbo out-of-plane dynamics during the first 
electrodynamic mission (it was extended to the out-of-plane oscillation 
damping issue and to an aocarate estimation of the tether bowing) is to be 
considered complete unless otherwise directed by NASA/M3PC, 

In the area of high resolul ?on slack tether studies, the SLACK2 
computer code will be soon available in its three dimensional formulation, 
Ve will start to simulate dynamics after the break in the most critical 
conditions (we look for inputs from NASA/MSFC), The study of the slack 
tether will be linked t s far as possible to the electroi3jnamic hazard 
investigation. 

In this area we plan to continue our analysis of the effects of the 
plasma contactor's failure, and to perform detailed calculations, using the 
SE7H00K code, of the voltages acquired, because of this failure, by each 
tether's end, with respect to the medium. 

In addition, we plan to analyze further the wire breakage cases. For 
these, we will construct an improved model of the truncated wire and with 
this new model we will replace the elongated prolate spheroid that we found 
to be inadequate. In addition, we will start the development of an 
approach for the evaluation of the effects of ionospheric plasma on the 
strength of the electric Meld in the vicinity of the truncated tip. 


u 
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The up-grade of the rotational dynamics model will be completed with 
the modelization of the major external perturbations, A now version of the 
plotting program in order to have outputs of more intuitive reading will be 
also worked out. Xf the debugging of the code does not beoome critical we 
should complete this level of upgrade by the end of the next reporting 
period. 

The development of safe operating area curves will be started as 
planned in the statement of work. We will appreciate any suggestions from 
NASA/MSFC related to what is considered first priority issue. 


